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1. Introduction 



Position coordinates and corresponding components of linear momentum of a free 
quantum particle are - - according to usual quantum mechanics - - complementary 
observables, a term coined by Niels Bohr in reference to measurements. The main 
expression of this property (for quantum mechanics on the real line) are Heisenberg's 
uncertainty relations AqAp > H/2, implying that there are no quantum states 
with arbitrarily narrow distributions of both conjugate variables q and p. If there 
were normalized states with Aq = 0, i.e. eigenstates of q, then Ap would diverge, 
hence measured values of p couldn't be predicted. However, due to continuity of 
spectra of q and p, the strict limits Aq = or Ap = cannot be physically 
attained. Complementarity thus means that quantum systems possess properties that 
are mutually exclusive: the observation of one of them precludes the observation of the 
other. A mathematical expression of complementarity of coordinate and momentum is 

\(p\q}\ 2 = const., (1) 

i.e. if we know everything about position, we know nothing of momentum, and vice 
versa. 

It is interesting and may have deep physical significance that pairs of complementary 
observables also exist for systems with finite-dimensional Hilbert spaces, as noted by 
J. Schwinger [T] (see also [21 E])- For such pairs precise knowledge of one of them 



implies that all possible outcomes of measuring the other one are equally probable, as 
exemplified in (pQ). 

If complementary observables have non-degenerate spectra, complementarity can be 
expressed in terms of the corresponding normalized eigenstates forming complementary 
bases. Then a measurement over one basis provides maximum uncertainty as to the 
outcome of a measurement in the other because all N possible outcomes will have equal 
probabilities 1/JV. The first attempt to use complementary bases in quantum state 
determination is due to Ivanovic [I] who provided explicit formulae for N + 1 such 
bases if N is an odd prime. The idea of using these bases for optimal quantum state 
determination was further developed by Wootters [5] and by Wootters and Fields [6] 
who called them mutually unbiased. In the latter paper they presented a construction of 
N + 1 mutually unbiased bases in an arbitrary prime power-dimensional Hilbert space 
and also demonstrated that they may serve as a maximal collection of measurements 
for optimal state determination. 

Our main concern are iV-level systems which provide basic models proposed 
for quantum information processing, since mutually unbiased bases find important 
applications there [7]. Their property that the outcome of a measurement in one 
selected basis gives no information about the possible results of measurements in all 
other mutually unbiased bases is of advantage for instance in key distribution protocols 
in quantum cryptography [8]. Observables with such a property for two-level systems 
(whose vectors are called qubits in quantum computation) are three Pauli matrices. 
Recently, <i-level quantum systems (with d > 2 and vectors called qudits) have come to 
a closer attention. It has been shown that such systems can be realized experimentally 
and quantum key distribution protocols using qudits have been introduced (see e.g. [9]). 
Since such protocols use mutually unbiased bases in dimensions higher than two, it is 
desirable to study constructions of mutually unbiased bases for higher dimensions, too. 
It is also known that complementary observables are useful in quantum state tomography 



2. Complementary observables and mutually unbiased bases 

Mutually unbiased bases in Hilbert spaces of finite dimensions are closely related to the 
quantal notion of complementarity. Let us start with two almost identical definitions. 

Definition 1 [3]. Two observables A and B of a quantum system with Hilbert space 
of finite dimension N are called complementary, if their eigenvalues are non- degenerate 
and any two normalized eigenvectors \ui) of A and \ vj) o/B satisfy 



Then in an eigenstate \ui) of A all eigenvalues b\, . . . , b^ of B are measured with 
equal probabilities, and vice versa. This means that exact knowledge of the measured 
value of A implies maximal uncertainty to any measured value of B. For the next 
definition note that the (non-degenerate) eigenvalues Oj of of A and bj of B are in fact 
irrelevant, since only the corresponding orthonormal bases \u{) and \vf) are involved. 

Definition 2 [5], [6] . Two orthonormal bases in an N -dimensional complex Hilbert 
space 




(2) 



{\ui)\i = 1,2,..., AT} and {\ Vj )\j — 1,2, ... , N} 
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are called mutually unbiased if inner products between all possible pairs of vectors taken 
from distinct bases have the same magnitude 1/y/N, 

\( Ui \ Vj )\ = -±= for all i,j E{1,2,...,N}. (3) 

In the above sense one may call two measurements to be mutually unbiased, if the 
bases composed of the eigenstates of their observables (with non-degenerate spectra) 
are mutually unbiased. Further, a set of bases is called mutually unbiased if every two 
different bases from the set are mutually unbiased. 

An important fact was proved in [6] (for further references see also [11]) that the 
upper limit to the maximal possible number of bases that can form a set of mutually 
unbiased bases in an ^-dimensional Hilbert space is N + 1: 

Theorem 1. In an N -dimensional Hilbert space, there cannot be more than N + 1 
mutually unbiased bases. 

Finally we remark that a criterion of equivalence of two pairs of mutually unbiased 
bases was formulated in [3]. 



3. Quantum structures in finite-dimensional Hilbert spaces 



Our starting point for quantum mechanics in the Hilbert space of finite dimension N 
is a model of quantum kinematics due to H. Weyl [12]. Its geometric interpretation as 
the simplest quantum kinematics on a finite discrete configuration space formed by a 
periodic chain of N points was elaborated by J. Schwinger [2]. In [T3J [T3] we proposed 
its group theoretical formulation based on Mackey's system of imprimitivity [13] which 
provides a group theoretical generalization of Heisenberg's commutation relations. For 
a recent review of the topic see [T6J [17] . 

In an iV-dimensional Hilbert space with orthonormal basis B = {\0),\l), . . .\N — 1)} 
we can establish a group generated by unitary operators Qn, Pn defined by the relations 

Q N \j)=<J N \j), j = 0,l,...,iV-l, (4) 
P N \j) = \j-l (mod iV)>; (5) 

here ujn is a primitive iV-th root of unity, e.g. ujn = exp(2iri/N). If B is the standard 
(or canonical) basis of C , the operators Pn an d Qn are represented by matrices 



Q N = diag (l,u N , 
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(7) 



In finite-dimensional quantum mechanics the unitary matrices Qn and Pn are analogues 
of exponentials of position and momentum in the continuous coordinate representation 
[12]. Namely, they fulfil an algebraic relation 

PnQn = ^nQnPn (8) 
which is analogous to the relation for Weyl's exponential form of Heisenberg's 



commutation relations. Further, P, 



N 
N 



Q 



Im, where Im is the N x N unit matrix. 
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This model has a simple geometric interpretation. The cyclic group Zjv = 
{0, 1, ... AT — 1} serves as a configuration space for iV-dimensional quantum mechanics. 
Elements of Z^r label the vectors of the basis B = {\0),\1), . . .\N — 1)} with the physical 
interpretation that \j) is the normalized eigenvector of position at j G Zjy. The natural 
transitive action of Z^r on Z N via addition modulo N is represented by unitary operators 
U(k) = Pjy. Their action on vectors \j) from basis B is given by 

U(k)\j)=P k N \j} = \j-k (modiV)) (9) 

The Zjv analogue of the Fourier transformation is the discrete Fourier 
transformation given by the unitary Sylvester matrix Sn with elements 

»=0W> = 4 (io) 



involving powers of u>n- The relations 

S N PnSn = Qni S n 1 QnSn = Pj^ 1 (11) 

show that the discrete Fourier transform diagonalizes the momentum operator, 
i.e. performs the transition from the coordinate representation to the momentum 
representation. 

The finite group generated by uj n , Q N and P/v 

ITv = WnQ'nIHI, i,j = 0,l,2,...,N-l} (12) 

consists of N 3 unitary matrices and is called the finite Heisenberg group [18] or the 
Pauli group [19]. It has been found useful in connection with mutually unbiased bases 
[2"T] . Note also that the set of N 2 unitary matrices [Q a N P b N \a, be {0,1, . . . , N - 1}} 
constitutes, as Schwinger [1] has shown, a complete operator basis of the Hilbert space 
of all complex matrices orthogonal in the sense of the inner product 

Tr ((Q^) t Q c N P d N ) = N8 ac 5 M for all a, b,c,d£ Z N . (13) 
4. Mutually unbiased bases for prime iV 

The question whether it is possible to attain the maximal number of iV + 1 mutually 
unbiased bases was answered in positive in [5J, but under a number theoretic proviso: 
a maximal collection of iV + 1 mutually unbiased bases exists in Hilbert spaces of 
dimensions equal to arbitrary powers of prime numbers. In this paper we will devote 
our attention to prime dimensions. 

In the particular case N = 2 one finds that the set of eigenvectors of the Pauli 
matrices forms a complete collection of three mutually unbiased bases: 

{ | 0> ,i J>} , m±ii>,MHi>\ m±M,wzjin (14) 



The construction of such a set of bases in higher dimensions can be understood as a 
generalization of this property. The formulae for iV + 1 mutually unbiased bases forming 
a maximal set for any odd prime dimension N were first given (not derived) in [1] ; we 
quote them according to [6j: 

\ v l 0) )j = 5 jk, 
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iV 



2-7ri l 

e * J . 



Here \vjf )j denotes the j-th component of the k-th vector in r-th basis, r = 0, 1, . . . , N. 
The first (r = 0) basis is the canonical basis, the last one (r = N) is its discrete Fourier 
transform. Mutual unbiasedness of the bases follows from the Gauss sums of number 
theory valid for p odd prime [20] 

p-i 

(15) 
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k=0 



e p 
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here a, b e N, a 7^ and a is not an integral multiple of p. 

A derivation of these N + 1 mutually unbiased bases has been given for any prime 
dimension N in [21J in terms of unitary operators Qn, Pn defined in Section 3: 

Theorem 2. Let N be a prime. Then the bases composed of eigenvectors of N + 1 
operators 



Qni Pn, PnQni PnQ 



2 

Ni 



PnQ 



AT-l 
N 



(16) 



are pairwise mutually unbiased and form therefore a maximal set of N + 1 mutually 
unbiased bases. 

In this paper we are going to give an alternative construction of a complete collection 
of iV + 1 mutually unbiased bases in a prime- dimensional Hilbert space by using the 
finite phase space related to the finite Heisenberg group. 



5. Finite phase space and its group of automorphisms 

In order to arrive at an independent proof of Theorem 2 we have to introduce the 
necessary group theoretical notions. 

First we need to establish a connection between the finite Heisenberg group and 
the finite phase space T N = Z N x Z N , N = 2,3, .. . rTSJ [TO]. The elements of 
the finite Heisenberg group were given in ([121) with l,i,j = 0, 1, . . . , N — 1. Its center 
Z(Un) is the set of those elements of Um which commute with all elements in Un, 

Z(U N ) = {(I, 0,0)\l = 0,1,..., N-l}. (17) 

Since the center is a normal subgroup, we can go over to the quotient group 11^ /Z(J1 N ). 
Its elements are the cosets labeled by pairs i, j = 0, 1, . . . , N — 1. The quotient 

group is then identified with the finite phase space T N = Z N xZ, N . To simplify notation, 
we shall denote the cosets corresponding to elements (i, j) of the phase space by Q % pi 
without subscripts N, 

Q l P 3 = {J N Q l N P 3 N \ 1 = 0,1,. ..,N-1}. (18) 

Note that all operators belonging to the same coset have the same eigenvectors because 
they differ only by multipliers u l N . 
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The definition of the phase space can be concentrated in the exact sequence [T5] 

i -> z{n N ) -> njv -> u N /z(n N ) -> = x z N -> 1 

expressing the relation of the projective representation of Z^ x Zjv to the central extension by Z(ILn) 
[12]. It is obvious that the correspondence 

: U N /Z{U N ) -> r w = Z w x Z w : Q«i* • • :>../). 
is an isomorphism of Abelian groups, since 

^ ((Q'i"') (V>')) = ((Q*^)) ((V'^')) = 
= + = (i + i'J+j')- 

We shall now focus on the group of automorphisms of the phase space T^. It was 
studied in |18j . However, we follow the approach of [22], where instead of cosets the one- 
dimensional grading subspaces of the Pauli graded Lie algebra gl(N, C) were considered 
and their transformations under the automorphisms of gl(N, C) were investigated. The 
subgroup of inner automorphisms of gl(N, C) was induced by the action 

tpx(A) = X~ l AX (19) 

of matrices X from GL(N, C). 

In the same vein we will concentrate on the automorphisms of the form ffl9|) . 
acting on elements of 1T7V, which induce permutations of cosets in U N /Z(I1 N ). Since 
the operators u^Q^P^ have the same spectra, the matrices X which induce the 
automorphisms ffl9|) are unitary. They can be understood as transformation matrices 
that transform a unitary operator (of the form u l N Q x P^) to a different basis, in which 
the operator is of the form uj m Q c P d . Thus X is a transformation matrix between two 
orthonormal bases. For explicit forms of matrices X see [18] (for N odd prime), but the 
results of [22] will suit better to our purpose. 

Automorphisms ip of the form f[T9~|) are equivalent if they define the same 
transformation of cosets in 11^ / Z{Yl^): 

^ Y r^^ x ^Y- 1 Q i P' J Y = X- 1 Q i P j X for all (i,j)eZ N xZ N . (20) 

Since the group Hn/Z(IIn) has only two generators — the cosets P and Q — condition 
ff20|) is equivalent to 

~^ x ^ Y- X PY = X~ l PX and Y^QY = X^QX. (21) 

If ipY induces a transformation of Un / Z (IIn) , then there must exist elements a, b,c,d G 
Zjv such that 

Y~ l QY = Q a P b and Y~ X FY = Q c P d . (22) 

It follows that to each equivalence class of automorphisms i/jy a quadruple (a, b, c, d) of 
elements in Zjy is assigned. We shall prove more, namely: 

Theorem 3. For N prime there is an isomorphism $ between the set of equivalence 
classes of automorphisms ipy o,nd the group SL(2, Z N ) of 2x2 matrices with determinant 
equal to 1 modulo N, 



$(if>y) = ( a c ^ , a,b,c,deZ N ; 



the action of these automorphisms on Un/Z(Un) is given by the right action of 
SX(2,Zjv) on the phase space = Zn x Zn, 

('',?) = (> J) fed)' < 23 > 
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Proof: To the composition of two automorphisms ip X) ipy corresponding to 
(a x , b x , c x , d x ) and ( CLy i by i Cy } dy), respectively, the product of matrices corresponding 
to ip x and tpY is assigned, as can be seen from 

{XYy 1 Q(XY) = {Y- 1 QY) ax {Y- l PY) bx = 

— Qiyax pbyax QCybx pdyb x _ Qa-xay+bxcy paxby+b x dy 

and similarly for P 

(XY)~ 1 P(XY) = Q CxaY+dxCY pcxby+dxdy 

Hence 

H^x^hr) = ^xM^Py) (24) 

and $ is an injective homomorphism. 

Now matrix elements a, b, c, d cannot be chosen arbitrarily. Consider the action of 

ifr- 

Y- 1 QY = Q a P b Y-'Q N Y = fiQ a N P b N , \fi\ = l, (25) 

Y~ X PY = Q c P d Y^PnY = XQ C N P^ |A| = 1. (26) 

By multiplying equation (1251) by equation ( 1261) once from the left and once from the 
right, we obtain 

P N Q N Y = fi\YQ c N P d Q a N P b N , Q N P N Y = ^\YQ a N P b N Q c N P d . (27) 
Using the commutation relation (|S]) we obtain 



CO 



- ad 



N 



fi\YQ% +c P b N +d = P N Q N Y = u^Q N P N Y = LU N l LU N bc ^\YQ a N +c P b N +d (28) 
leading to the condition 

oo N ad = uJ N bc -\ (29) 
It will be fulfilled if and only if ad — be = 1 (mod N), i.e. 

det (* b A=l (mod N). (30) 



c d 

This means that to every ipy acting on Un/Z(Un) a matrix from SL(2, Zjv) 

a b 



c d 



is assigned. Now to every coset from njv/Z(II/v) an element of the phase space 
Zjv x Zjv was associated. So finally we check that the action of ipy on Q l Pi is given by 

(i',/) = Q*'^' = i) Y (Q l PJ) = Y~ 1 Q i P j Y = 

= Y- l Q l YY- l P j Y = Q^ocpib+jd = ( ia + j c ^ b + j d ^ 

and this means that the transformation of can be written as the right action of 
SL(2,Z N ) on Ztv x Z n 

a J) = (u) ( ° J ) ■ (3D 

Finally observe that mapping $ is an isomorphism, since 
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follows from 



(i'j') = (ij) if and only if X^Q'P^X = I^Q^Hn. □ 
We conclude this section with 

Lemma 1 [23J. The right action of SL(2, Zn) on the phase space Zn x Zn does 
not change the determinant of a matrix composed of components of two vectors from 

Zj\r X Zjy. 

Proof. Consider two vectors (i, j) and (k, I) from ZnxZn and the matrix 

Since the action of A is given by (i',f) = (i,j)A and (k',l') = (k,l)A, and because 
det A = 1, one immediately gets the result 

det(£ £)=d*((j[ f )>) = 

= det(* j)de tj 4 = det(° j). □ (32) 

Remarks. For A" prime the right action of SL(2, Zn) on the phase space Zjv x Zj\r 
has exactly two orbits — the single point {(0, 0)} and On = Zn x Ztv\{(0, 0)} consisting 
of A^ 2 — 1 points. The stationary subgroup of the point (1,0) from On is the Abelian 

subgroup { ( ^ J J \b — 0, 1, . . . , N — 1} of order N. Hence the order of SL(2, Zn) is 

N(N 2 — 1). Further, according to Lemma 1, the determinant fl32l) is an invariant of the 
right action of SL(2,Zn) on r^r X ^n- Let us note that SL(2,Z^r) transformations of 
the finite phase space were also studied in [24|. 



6. New construction of the maximal set of mutually unbiased bases for A^ 
prime 

In this section the finite phase space and its transformations of the form (Tl9|) will be 
used to introduce an interesting algebraic structure that proves the existence of A^ + 1 
mutually unbiased bases for prime N, thus providing an alternative approach to their 
construction. We shall exploit the fact that, for prime N, Zn is a finite field, i.e., there 
is also a multiplicative group structure modulo A^ in Z* N = Zjv\{0}. 

Our construction starts with the partition of the finite phase space r^v = Zn x Zn 
into equivalence classes [(i, j)} defined by the equivalence relation ~ if 

there exists r e Z* N such that (i',j f ) = (ri,rj), where the multiplication is understood 
modulo N. 

We exclude the trivial class [(0,0)] containing only (0,0). Then the orbit 
On = Zn x Ztv\{(0,0)} is decomposed into N + 1 classes [(1,0)] and [(«, 1)] where 
i = 0,1, . . . , N — 1. The fact that Zn is a field for prime A^ is crucial in the proof 
that every element of the orbit On belongs to some class. Since each class has A^ — 1 
elements, this decomposition contains A^ 2 — 1 elements in total, with the only element 
(0, 0) not included. 

If an element is of the form (0, i) or (i, 0), then it is obvious that it belongs to classes 
[(0,1)] or [(1,0)], respectively. An element of the form i,j = 1,2,..., AT — 1, 

will belong to the class [{k, 1)] where k G {1,2, . . . , N — 1} is the solution of kj = i 
(mod N). The existence and uniqueness of such k is guaranteed by the fact that Zn is 
a field for A^ prime. The partition into classes can be visualized in the following table 
of all elements in Z N X Zjv\{(0, 0)}: 
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(0,1) 


(0,2) 


■ ■ (0,N-1) 


1 
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(1,0) 




(1,1) 


(1,2) • 


(~\ 1\T 1 ^ 
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(2,0) 


(2,1) 


(2,2) • 


■■ (2,N-1) 


N-l 


(N-1,0) 


(N-1,1) 


(N-1,2) • 


• • (N-1,N-1) 



Here every element corresponds to a coset Q l PK All operators in the same 
coset differ just by a complex multiplier. Every multiple (ri,rj) of a vector by 
r G {1, 2, . . . , N — 1} will therefore correspond to the coset Q"P rj . Because of relation 
(EJ) it is obvious that operators (<5 i -P J ) r and (Q n P r i) belong to the same coset. An 
important consequence is that elements (ri,rj), r = 1,2, . . . , N — 1, correspond to 
commuting operators, hence have the same eigenvectors. Thus we have proved 

Lemma 2. If N is a prime, then there are exactly N + 1 classes of elements from 
On = x 2W\{(0,0)}, each class containing N — 1 elements. All elements of the 
same class correspond to commuting operators with the same eigenvectors. 

We will now demonstrate that the bases composed of the eigenvectors of two 
different operators corresponding to elements from distinct classes in On are mutually 
unbiased. 

Theorem 4. Let N be a prime and let (a, b) and (c, d) be two elements from 
On = %>n x ^A r \{(0,0)} which belong to distinct classes [(a, b)) ^ [(c, d)]. Then the 
bases composed of eigenvectors of the operators from the corresponding cosets Q a P b and 
Qcpd are mu tually unbiased. 

Proof. The first step is to show that the bases composed of eigenvectors of Qn and 
Pn are mutually unbiased. This follows directly from equation ( JTTj) . Namely, 

PnS n \j) = S n Qn\j) = uJ N S N \j), (33) 

where \j) is an eigenvector of Qn, so Sn\J) is an eigenvector of Pn- Further, because of 
f fTUj) . the inner product of \j) and S^k) has absolute value 

\(\ 3 },SN\k))\ = \(j\S N \k}\ = 

Hence if we have two elements where one belongs to the class [(1,0)] and the other 
to the class [(0, 1)], then we already know that their corresponding bases are mutually 
unbiased, because they are composed of eigenvectors of Qn and P N , respectively. 

Because of the partition of On it is sufficient to consider now only the case of two 
distinct elements (a, 1) and (b, 1), with a, b e {1, 2, . . . , iV — 1}, a ^ b. We are going to 
show that the bases of eigenvectors of the corresponding operators Q n Pn and Q n Pn are 
mutually unbiased (hence also the eigenvectors of powers of these operators). According 
to Theorem 3, to unitary operators X that permute the cosets in the Heisenberg group 

X^QipiX = Q l 'P j ' 

matrices from SL(2, Z^r) are assigned. Conversely, to every matrix from SL(2, Z^r) there 
is an equivalence class of unitary operators which induce the same permutation of the 
cosets. In this sense a special unitary representation of SL(2, Z N ) was described in [18J. 



UJ 



N 



(34) 
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We will now show, if a ^ b, then there exists a matrix A from SL(2, Z^r) such that 

(a,l)A=(5,0) and (b, 1)A = (0, 6) 

If we indeed can find such a matrix, then there exists a corresponding unitary operator 
X such that 

X~ 1 Q a PX = Q a , X _1 Q b PX = P S , 

hence the eigenvectors of Q a P and Q b P can be expressed as X\j) and XSat|&), 
respectively. According to (|34|) modulus of their inner product is 

|(Xb-),^|*»| = |(b->,fiW|*»|= ' 



'JV 

proving that these bases are mutually unbiased, too. 

To prove the existence of a unique matrix A £ SL(2, Zat) with the desired properties 
we apply Lemma 1 implying 

det ( Q b \ J = a - b (mod N) = det ~ J = ab (mod AT) 

and we select a, b £ Z^ such that ab = a — b (mod AT). Equivalently we look for a 
matrix C = A^ 1 £ SL(2,Ztv) producing the inverse transformation 

(a,0)C = (5,0)(^ £)=(a,l), 

(0,&)<7=(0,S)(° = (M). 

This gives us the following equations to compute the elements of C: 

a(3 = 1 (mod N), (35) 

aa = a (mod N), (36) 

67 = 6 (mod A/"), (37) 

65 = 1 (mod AT). (38) 

The fact that AT is a prime guarantees that each of these equations has unique solution 

a (3 



in 7h N . Having the values of entries a, ft, 7, 8, we still need to check that C 



7 8 



belongs to SL(2,Zat). By multiplying equations fl36l) and (1381) and subtracting the 
product of (13"7j) and (1331) we obtain 

ab(a8 — (3j) = a — b (mod AT). 

Since ab = a — b (mod AT) we have 

det C = a8 - fa = 1 (mod N), 

verifying that C indeed belongs to SL(2,Zjv)- The inverse matrix A = C^ 1 = 

_^ a ) wm then transform pairs (a, 1) and (6, 1) into (a, 0) and (b, 0), respectively: 

(a,!)^ 1 = (S,0), (M)C _1 = (0,6). 
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To complete the proof, it is easy to see that for pairs (6,1), (1,0) and (6,1), (0,1), 
6 = 1, 2, . . . , N — 1, there exist unique transformation matrices from SL(2,Ztv) such 
that 

(6,1)4(6) = (0,1), (1,0)4(6) = (1,0) 

and 

(6,1)4(6) = (6,0), (0,1)4,(6) = (0,1) 

Hence pairs of bases composed of eigenvectors of pairs of operators Q b P, Q and Q b P, P 
are mutually unbiased. Thus we have shown that there exist N + 1 mutually unbiased 
bases in a Hilbert space of prime dimension N. We have therefore reached the same 
conclusion as [21] . In our case the mutually unbiased bases are composed of eigenvectors 
of operators 

Qn, Pn, QnPn, Qn^n, ■ ■ ■ , Qn 1 Pn, (39) 

while the operators in Theorem 2 are only modified using (jSJ). 

Note that one could make a different choice of representatives of the classes forming 
the partition: we could have alternatively used e.g. the pairs (1, a), a = 1, 2, . . . , N — 1, 
instead of (a, 1), and the mutually unbiased bases would be given by bases composed of 
eigenvectors of operators 

QniPniQnPn-iQnPni ■ ■ ■ iQnPn l - 

To provide a constructive proof we should give an explicit way to construct the 
bases out of the canonical basis £>. Let us denote the bases composed of eigenvectors of 
(El by 

# = #(1,0), #(0,1)5 #(1,1)5 #(2,1), •••) #(7V-1,1)- 

We know that the map #(i,o) — * #(o,i) is implemented by the unitary operator Sn- 
The next step #(o,i) —* #(i,i), leaving #(1,0) intact, clearly corresponds to the above 
transformation matrix 4( — 1) from SL(2, Zjv). Its iterations will generate further steps. 
A unitary transformation Dn which implements 4( — 1), 

D^QnP>n = Qn, Dn PnDn = £n QnPn, 
can be taken from [22J; it is diagonal, 

D N = diag (d , d 1: . . . , djv-i), dj = e N J uj^, 

where En = 1 if iV is odd, En = ^Ju^N if iV is even. In this way we arrive at a sequence 
of unitary maps 

10 Sn yj Dn Djv y2 ®n £>JV jo 

"(1,0) — > £>(0,1) —> "(1,1) —> "(2,1) 0(jv-l,l), 

and the composite unitary operators D b N SN, 6 = 0, 1, . . . , N — 1 will produce all the 
bases starting from the canonical one. □ 

Example N = 2. The phase space T 2 consists of 4 elements (0,0), (1,0), (0,1), 
(1, 1). The group SL(2,Z 2 ) (also known as the group of invertible 2x2 matrices over 
the simplest finite field F2) with 6 elements 

(I 0\ (0 1\ fl 0\ (1 l\ (1 1\ (0 l\ 

\o 1 O/V 1 1 /'V° 1 / ' V 1 °/'V 1 1 / ' 



m» = ( i ! ) 

M0) = ( J V ) 
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acts transitively on the orbit {(1, 0), (0, 1), (1, 1)}. Unitary operators transforming the 
bases flH} 

#(1,0) £(o,i) 

are 



7. Concluding remarks 

The question whether the maximal number N + 1 can be attained for a composite 
dimension N, where N is not prime nor a power of a prime, still remains an open 
problem. The answer is not known even for the simplest case N = 6 where N + 1 = 7. 
A simple argument [21] leads to 3 such bases. Some numerical attempts to find further 
mutually unbiased bases were not successful. So it remains unclear whether it is indeed 
possible to reach the maximal number N + 1 of them for the Hilbert space of dimension 
6 and in other composite dimensions as well. 

Although the relation between the eigenvectors of Q % N P 3 N and mutually unbiased 
bases was observed e.g. in [2T], the relation between the decomposition of the phase 
space Ztv x Ztv whose elements correspond to cosets in the finite Heisenberg group 
and the existence of mutually unbiased bases has been left unnoticed so far. We were 
thus able to give an independent constructive proof of Theorem 2 using group theory. 
However, the proof heavily depends on properties that are a consequence of N being a 
prime. The extension of our proof to the case of Galois fields with prime powers N will 
be subject of a future publication. It might also be interesting to investigate whether 
our procedure would provide better insight in the problem of existence or non-existence 
of the maximal number of mutually unbiased bases in composite dimensions. 
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